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Abstract: The dynamics of an optically trapped particle are often 
determined by measuring intensity shifts of the back- scattered light from 
the particle using position sensitive detectors. We present a technique which 
measures the phase of the back- scattered light using balanced detection in 
an external Mach-Zender interferometer scheme where we separate out and 
beat the scattered light from the bead and that from the top surface of our 
trapping chamber. The technique has improved axial motion resolution over 
intensity-based detection, and can also be used to measure lateral motion 
of the trapped particle. In addition, we are able to track the Brownian 
motion of trapped 1 and 3 fim diameter beads from the phase jitter and 
show that, similar to intensity-based measurements, phase measurements 
can also be used to simultaneously determine displacements of the trapped 
bead as well as the spring constant of the trap. For lateral displacements, 
we have matched our experimental results with a simulation of the overall 
phase contour of the back- scattered light for lateral displacements by using 
plane wave decomposition in conjunction with Mie scattering theory. The 
position resolution is limited by path drifts of the interferometer which we 
have presently reduced to obtain a displacement resolution of around 2 nm 
for 1.1 jim diameter probes by locking the interferometer to a frequency 
stabilized diode laser. 
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1. Introduction 

A single microparticle trapped by optical tweezers can be used as a micro-probe which can 
be used in diverse applications. In some of these applications, the object of interest (such as 
DNA/RNA, molecular motors such as myosin and kinesin, single bacteria, etc.) is attached to 
the micro-probe which is held controllably in a single trap (or in some cases multiple probes 
trapped separately), and its motion carefully monitored to yield information about the dynamics 
of the object of interest including micro-forces and torques [Q]EIEl|4l[5]|6]], while in other cases, 
the probe itself is used to reveal interesting information about surface topographies with nm 
precision (7), Brownian motion (HO, and to study fundamental statistical physics phenomena 
ifTOl . The probe motion is manifested in changes in the forward and backward scattering pat- 
terns of a detection laser incident on the trapped probe - the scattering patterns being typically 
imaged on position sensitive detectors (PSD) or quadrant photodiodes (QPD). Backscattered 
detection is often preferred over forward scattering due to its relative independence of the mor- 
phology of the trapping apparatus, and has been shown to resolve displacements even in the 
pm regime ifTTTl . However, the high sensitivity obtained in displacement sensing using PSDs or 
QPDs is for probe motion in the radial direction, i.e. the in the direction perpendicular to the 
trapping beam. For motion in the axial direction, the resolution is much lower since in this case, 
the intrinsic resolution of QPD or PSD is no longer obtained, the motion being determined by 
measuring the change in the total amount of light incident on all the four quadrants of the QPD 
(or all the pixels of a PSD), and not the difference between pairs of quadrants (pixels for a PSD). 

In this paper, we experimentally demonstrate a technique where the phase of the back- 
reflected beam is measured as a function of probe motion instead of the intensity. It has been 
shown earlier lfT2l that phase measurement of the interference pattern formed by superposition 
of the backscattered light from the probe and the glass trapping chamber can be used to quantify 
probe axial motion. We improve the sensitivity of axial motion detection by improving the sig- 
nal to noise of the interference signal by using balanced detection in an external Mach-Zender 
interferometer configuration. In addition, we show that even the radial motion of the probe 
causes a change in the phase of the interference signal, which can thus be used to calibrate 



radial motion as well. To calibrate the change of phase with radial motion, we perform a sim- 
ulation using plane wave decomposition in conjunction with general Mie scattering theory to 
find out the phase shift as a function of probe displacement in the radial direction. Our method 
is also sensitive to Brownian motion of the trapped probe - which we are able to quantify for 
a given trap stiffness from the phase jitter of the interference pattern. Also, as is the case in 
intensity -based detection, a fourier transform of the phase jitter yields the corner frequency of 
the probe, and thus the spring constant or stiffness of the optical trap fT3l . 

2. Materials and methods 

2.1. Experimental apparatus and methods 

Our optical tweezers apparatus is developed around a Zeiss Axiovert Al inverted fluorescence 
microscope as shown in Fig. [I] A 100X, 1.4 N.A. oil immersion microscope objective (Zeiss 
plano-apochromat, infinity corrected) is used to couple the beams into the sample chamber. 
Polystyrene beads of size 1.1 or 3 jim (which serve as our micro-probes) immersed in wa- 
ter are trapped using a focused single transverse mode 1064 nm solid state laser (Lasever 
LSR1064ML) beam of 600 mW power and specified M 2 quality factor of 1.2. An acousto- 
optic modulator (AOM) is inserted in the back focal plane of the trap to move the position of 
the probe radially in a controlled manner. A pair of convex lenses are used to image the plane 
of the AOM on to the back aperture of the microscope objective to ensure that any angular 
deflection at the AOM gets directly mapped to an angular deflection at the back aperture of the 
objective without the beam walking off. The efficiency of the AOM is around 50% for the first 
order, so that the maximum power available for trapping is around 180 mW after other losses in 
the coupling optics. Detection is performed using a separate laser beam (Toptica DL 100, wave- 
length 780 nm, max power 75 mW) whose power is kept low (around 1 mW finally) enough 
so as to not modify the trapping potential. The trapping and detection beams are combined at 
the input of the microscope using a polarizing beam splitter. The detection laser beam is back- 
reflected from the probe and collected at the microscope side port using a dichroic mirror which 
has very low transmission (but very high reflection) at 1064 nm but about 50% transmission at 
780 nm. Typically, this back-reflected light is incident on a miniature quadrant photo detector 
ifTH to quantify the radial motion of the probe, and determine the power spectrum of the probe 
motion for trap calibration. The sample chamber consists of a glass slide coated with gold (top 
surface of the chamber) and a cover slip of 160 jim thickness (bottom surface facing the ob- 
jective), with the polystyrene probe suspension sandwiched between. The gold coated slide is 
used to enhance the back-reflection from the surface which is required to compensate for the 
low reflectance of glass given that the power of the detection laser is quite low. About 25 fll of 
the sample is used at a dilution of 1:10000. 

The light at the microscope back-focal plane consists of back- scattered light (at 780 nm) 
from the trapped probe and unscattered light from different regions of the sample chamber. 
However, in our experiments, we trap the bead close to the top slide and separate out the scat- 
tered and unscattered components by a confocal arrangement consisting of a spatial filter used 
in combination with two apertures (Fig. [I]). The spatial filter comprises of a focusing lens and 
a pinhole of diameter 10 /im. The focal length of the lens is carefully chosen so that the waist 
size at the focal spot is close to the pinhole diameter. The lens focuses the backscattered in- 
tensity pattern on the pinhole to produce an Airy pattern at the pinhole output. The diameter 
of the detection laser is also kept slightly larger than the probe so that the edges of the beam 
pass by the probe and are reflected by the top slide directly. The central region of the Airy pat- 
tern thus corresponds to light backscattered from the trapped probe (which is close to the focus 
of the detection laser), while the comparatively diffuse outer rings correspond to the reflected 
light from the gold surface. Then, using two more apertures as shown in Fig. [T] a portion of 
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Fig. 1: Schematic of the experiment. Key: PD1: Photodiode 1, PD2: Photodiode 2, BS1 : Beam 
splitter 1, BS2: Beam splitter 2, Signal Gen.: Signal Generator. The pinhole generates an Airy 
pattern - the central portion of which contains scattered light from the trapped probe, while 
the diffuse ring contains unscattered reflection from the top gold slide of the sample chamber. 
Aperture 1 picks off the probe signal while Aperture 2 picks off a portion from the ring. The 
trapping laser is at 1064 nm and the detection laser is at 780 nm. 

the light from the outer rings can be separated from the central pattern. We proceed to beat the 
two separated components in a Mach-Zender interferometer where the path length of one of the 
arms is modulated by a piezo mirror to obtain interference fringes on both photodiodes PD1 
and PD2 (Thorlabs DET110). The signals on the photodiodes are out of phase by 180°, so the 
difference gives twice the individual signal amplitude. This is well known as balanced detection 
and improves our fringe contrast, and thus the sensitivity of the phase measurement by a factor 
of two over standard backscattered interferometry. The difference signal is obtained using a 
lock-in amplifier (SR630), which also measures the phase of the output signal with reference 
to the phase of the driving signal to the piezo mirror. The final output of the lock-in is shown in 
Fig.|2ja). Now, the phase of the output signal from the interferometer changes when the probe 
moves with respect to the top slide, with a 2n phase change for axial probe displacement of a 
unit wavelength of the detection laser. However, the estimation of phase change for axial mo- 
tion of the trapped probe with respect to the top slide needed a calibration of the axial distance 
moved. This could not be performed in our apparatus due to the lack of a three dimensional 
piezoelectric motional stage for our microscope. The probe was instead moved transversely 
(i.e. in the radial direction) by the AOM, and we used the fact that there is a phase change 
even in this case since the axial depth actually varies due to the curvature of the bead. Fringes 
obtained in the interferometer due to transverse motion of the bead are shown in Fig.[2jb). Ex- 



periments were performed with probes of diameter 1 . 1 and 3 fim. For the distance calibration, 
we used the pixel-to-distance calibration tool in our microscope camera (Axiovision) software 
and used the diameter of the probes as reference. Phase measurements as a change of probe 
position are shown in Fig. [5] Each data point was calculated by averaging over 10 independent 
measurements. 
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Fig. 2: (a). Balanced detection signal obtained by subtraction of the two out-of-phase signals 
from PDl and PD2, the photodiodes kept in the two arms of the Mach-Zender interferometer, 
(b) Typical fringes obtained in the phase measurement as the AOM voltage is changed to dis- 
place the trap transversely and cause a phase shift in the backscattered signal that is captured in 
the balanced detection output. 



2.2. Theoretical simulation 

The phase change of the backscattered signal due to axial motion is quite straight-forward 
(units of 27T per half wavelength change of axial distance). On the other hand, the same change 
when the probe moves radially is non-trivial. A theoretical estimate was thus required to match 
the change in the phase of the back- scattered pattern with experimental results. While there 
exists literature on the intensity of the backscattered light field in optical tweezers (HJ, we are 
not aware of any study of the phase profile of the backscattered pattern. Furthermore, even 
Ref. fT5l does not take into account the effect of the back-reflection from the top slide of a 
trapping chamber in calculating the intensity pattern. This is even more critical in determination 
of the phase profile, since the back-reflected light would interfere with the light scattered from 
the bead itself, thus modifying the phase profile of the overall backscattered signal. Therefore, 
a model was required to enable the understanding of the alterations in phase contour as the 
trapped microsphere was moved transversely. To develop the model, we used a variant of the 
Angular Spectrum Method (also referred to as vectorial Debye diffraction theory or Debye 
integral) lfT6l to calculate the electric field distribution. 

In our model, we considered a x-polarized Gaussian beam of light incident on the micro- 
sphere, having the form 

E(x,y,z) =E W ° exp ( ^ — exp ( -ikz - ik— r + i£(z-z w )\ i (1) 

w{z-z w ) \w z {z-z w )J \ 2R(z-z w ) J 

where r = \/x 2 +y 2 is the radial distance, Eo is the peak intensity which was set to unity, wo is 
the size of the waist or the tightest spot of the beam, z w is the position of the waist relative to the 
z = plane, w(z) is the size of the waist of the beam in the specified z plane, with w(0) = wo, 
R(z) is the radius of curvature of the phase front in the specified z plane, and is the Gouy 
phase shift. 

As is well known, a Gaussian beam can be decomposed into an infinite number of plane 
waves with appropriate weight factors. Each of these plane waves interact with the particle in 



accordance with the theory of Mie scattering (we consider the case when the wavelength of 
the light is comparable with the size of the scatterer). This decomposition was performed by 
applying a two dimensional fast fourier transform on the incident light field at the position of 
the microsphere. 
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Fig. 3: The theoretical model 
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Here, M and N are the sizes of the 2-dimensional array for the image of the cross section of 
the incident Gaussian beam. In each of the k- vector directions, a plane wave of unit magnitude 
was assumed to be generated which eventually scattered off the microsphere. The plane wave 
was assumed to be x polarized in it's own frame of reference (defined as x'). 

Ei = e ikz e[ (3) 

For such an incident plane wave, the scattered wave is of the form 
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and a m b n are coefficients of scattering, while N^ J n and M^ n are vector spherical harmonics 
with m=l |[T5lh given by the following expressions 

N fin =cos0 w(n+l)sinfl?fo(cosfl) — ^ e r 

P 

-f cos (j)T n (cos G) eg 

m^ [ph " ){p)] 
— sm(j) 7T„(cos 6)hn — (6) 

P 

and 

M oil = cos0 7t n (cos0)h£\p)e e - sin0T w (cos0)4 1) (p)e^, (7) 
In these expressions, 



P =kr, 



with being the spherical Hankel function. The scattering coefficients were calculated from 
the boundary conditions, and are given by the following when the permeability of the medium 
and microsphere are the same. 



where, 
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Here ft m is refractive index of the medium, n p is the refractive index of the microsphere and \jf n9 
£ n are the Ricati- Bessel functions. 

The scatter patterns corresponding to each plane wave were transformed from r, 6 and (j) basis 
to the x f , y' and z! basis of the plane wave frame using the following transformation matrix 



(10) 



The components of the scattered field in the x f , y' and t! basis were then transformed into the 
x, y and z basis of the lab frame using an appropriate coordinate transformation, and eventually, 
scattered contributions from all the plane waves were added up to form the final backscattered 
pattern as shown in Eqn.fTT) 
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Once the back- scattered field was calculated at a certain location, the field of the light re- 
flected from the top glass slide was also estimated at the same location by using Eq.[T]and the 
difference in the two fields calculated. The phase contour of the resultant field was then taken. 
In order to simulate the shift in the longitudinal position of the microsphere from the waist of 
the incident detection laser beam, we used an aperture function that was moved similarly across 
the microsphere. The size of the aperture function was taken as the size of the microsphere it- 
self. As the aperture function was moved, the FFT in the reciprocal plane shifted as well. The 
greater the offset in the radial direction, the more the shift of the FFT from the center in the 
reciprocal plane as demonstrated in Fig [4] Therefore, this also produced a scatter pattern which 
was shifted from the center in the transverse direction. Next, the change in phase of the resul- 
tant field (superposition of scattered and reflected) was determined as a function of aperture 
offset by comparing the phases of the fringes at different offsets. Results of the simulation are 
provided in Fig [5] 



Fig. 4: FFT of the incident Gaussian as the aperture function is offset in the x direction by (a) 
-1 jUm (b) jim (c) 1 jUm. The FFT plane is assumed to be 2 jim from the focus in z direction. 
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Fig. 5: Experimental and simulation data for phase change for light scattered off 1.1 jUm and 
3/im microspheres (probes) for known travel in the radial direction. In the experiment, the 
probes are moved in the radial direction by the AOM, while in simulation, we translate the 
aperture across the probe. The error bars in the experimental data signify la standard devi- 
ation. The standard deviation is predominantly due to drifts in the interferometer path length, 
but at low trapping powers, the Brownian motion of the trapped probe also contributes. The 
data point for the highest displacement of the 3 jim probe has a large error bar since the back- 
scattered signal to noise itself was low with the bead having been displaced significantly from 
the detection laser. 



3. Results 

A comparison of simulation and experimental results are shown in Fig. [5] The experimental 
data points have error bars that are due to: a) path drifts of the interferometer, and b) Brownian 
motion of the trapped probes. The path drift is a limiting factor of our experimental technique 
in the sense that it determines the minimum displacement resolution achievable. In the data 
shown here, the path drifts were of the order of 15 degrees that corresponded to around 75 nm 
in terms of displacement. Note that the transfer function between displacement and phase shift 
was determined by a straight line fit to the theoretical data, the slope of which was 4.9 ± 0.8 
/im/deg for 1.1 /im diameter probes. However, the path drifts can be reduced by locking the 
interferometer to a frequency stabilized diode laser by a simple side of fringe locking technique. 
By this, we managed to reduce the phase drifts to around 40 mdeg over averaging times of 100 
ms, so as to give a possible displacement resolution of about 2 nm at a bandwidth of 10 Hz for 
1.1 fim diameter probes. The phase measurement technique can then be implemented by either 
scanning the interferometer by a master-slave technique 1 17 ], or by using a second frequency 
stabilized laser that can be scanned independently to produce fringes in the interferometer. It 
is interesting to note that the radial displacement resolution can be higher for smaller probes 



since the curvature of such probes would be higher resulting in greater phase change for smaller 
displacements. However, the backscattering cross section would also be lower in these cases, 
resulting in reduced signal to noise that would serve as a check to the resolution achievable. 
The interferometer path drifts could be reduced by using temperature stabilized or fiber-based 
cavities that could thus increase the resolution of phase measurement and even achieve sub-nm 
resolution in displacement sensing in axial as well as radial probe motion. 

The Brownian motion of the trapped probes was also detected in our measurements. This 
was evident from the fact that the standard deviation in the data for optical powers of 47 and 
120 mW comes out to be 19.0(0.3) and 16.1(0.2) deg respectively, the uncertainties being the 
statistical uncertainties after averaging 10 times. While the standard deviation due to path drifts 
is 15(0.2) deg, the enhanced standard deviation would have to be due to another independent 
process, which from the inverse dependence on power, could only be Brownian motion. From 
the standard theory of errors, we could thus use the fact that a} ot = oj ri f t + cr^ m , where of ot is 
the total measured variance, crj^and <j£ m being the variances due to path drifts and Brownian 
motion respectively. Since errand are known from experiment, cr^ m , and thus cr^ m - the 
standard deviation due to Brownian motion can be determined. This comes out to be 55(5) 
nm and 25(5) nm for optical powers of 47 and 120 mW respectively. These results could now 
be compared to theoretical estimates of Brownian motion at the given trapping powers. From 
Ref. fT8lL the amplitude of Brownian motion of a trapped probe of radius r is given by the 
expression 



Here, K is the spring constant of the optical trap, kg is the Boltzmann constant, T is the tem- 
perature, r] is the dynamic viscosity of water, and t ave is the averaging time. The spring constant 
can be determined by a measurement of the corner frequency f c of the trapped particle from 
the expression 



where, yo = 6nrr[ . A standard measurement of corner frequency can be performed by the power 
spectrum method fT3ll . It is well known that a trapped probe executing Brownian motion obeys 
a simplified Langevin equation so that the power spectrum of the probe motion is a Lorentzian. 
Now, a position sensitive detector or a quadrant photodiode is used to record the displacement 
of the probe using a detection laser so that the power spectrum can then be obtained. In our 
technique, the phase has a linear relationship with the displacement - thus, a fourier transform 
of the jitter in the phase yields a similar power spectrum as is shown in Fig. [6] A Lorentzian 
fit to the data yields a corner frequency of around 25 Hz, which leads to a trap stiffness of 
1.3 pN//im for an optical power of 47 mW. Then, using a = 0.55/im, an integration time of 
t ave of 8.3 ms, temperature T of 300K, viscosity 7] of 0.008 Poise, one obtains the theoretical 
estimate of Brownian motion as 49(2) nm, the error being due to the fit error in the corner 
frequency. This is within la of our experimental estimate of Brownian motion at the same 
power. A similar measurement at an optical power of 60 mW yields a theoretical estimate of 
21(2) nm, which is again within lcr of the experimental measurement. At higher powers, the 
path drifts dominated the phase jitter and it was not possible to separate out the Brownian 
motion of the probe. As a check of the corner frequency measurements using the phase jitter, 
we performed power spectrum measurements of the displacement jitter of the probe at the same 
power levels using our quadrant photodiode detection system |[T4lL and found an agreement to 
within 10%. We can therefore conclude that our phase measurements can be used to determine 
probe displacements to the level of a few nm, as well as trap stiffness for an optically trapped 
probe in photonic force microscopy. The technique is more sensitive to axial displacement, 
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f c = K/(27iy ) 
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however, we can use it to determine radial motion as well with the help of a transfer function 
determined from experiment and verified by simulation. 
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Fig. 6: Power spectrum obtained by fourier transform of the phase jitter for a trapped probe of 
diameter 1 . 1 jim for an optical power of 47 mW. The data is fitted to a Lorentzian using IGOR 
fitting software. The fit parameters are shown, with the corner frequency given by f c = y/(B) 
coming out to be about 25 Hz. 



4. Conclusion 

In conclusion, we have developed a phase sensitive interferometric technique for simultane- 
ously measuring probe displacements and trap parameters in a photonic force microscopy setup. 
The technique is based on balanced detection of the output of two arms of a Mach-Zender in- 
terferometer set up using the backscattered light from the trapped probe and its environment, 
and reduces the inherent lack of sensitivity in standard back focal plane interferometry for axial 
probe displacements. We have extended the technique to measure radial probe displacements 
and matched our experimental results for radial motion with a theoretical simulation. The sim- 
ulation was performed using plane wave decomposition in conjunction with Mie scattering 
theory to find out the phase distribution of the backscattered signal due to radial probe motion, 
also taking into account the effects of a sample chamber. In addition, our technique is sensitive 
to Brownian motion, and can be used as such in any experiment to determine Brownian motion 
optically. The displacement sensitivity is limited mostly by the path drifts of the interferometer, 
that we have controlled presently to a level where we can achieve a resolution of around 2 nm 
for 1.1 fim diameter probes. This could be improved further by using temperature stabilized or 
fiber-based cavities, so that the capabilities of the technique can be extended to achieve sub-nm 
resolution for probe displacement in photonic force microscopy. 
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